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We consider stationary rotating black holes in SU(2) Einstein- Yang-Mills theory, coupled to a 
dilaton. The black holes possess non-trivial non-Abelian electric and magnetic fields outside their 
regular event horizon. While generic solutions carry no non-Abelian magnetic charge, but non- 
' Abelian electric charge, the presence of the dilaton field allows also for rotating solutions with no 

non-Abelian charge at all. As a consequence, these special solutions do not exhibit the generic 
asymptotic non-integer power fall-off of the non-Abelian gauge field functions. The rotating black 
hole solutions form sequences, characterized by the winding number n and the node number k of their 
gauge field functions, tending to embedded Abelian black holes. The stationary non-Abelian black 
hole solutions satisfy a mass formula, similar to the Smarr formula, where the dilaton charge enters 
instead of the magnetic charge. Introducing a topological charge, we conjecture, that black hole 
solutions in SU(2) Einstein- Yang-Mills-dilaton theory are uniquely characterized by their mass, their 
angular momentum, their dilaton charge, their non-Abelian electric charge, and their topological 
charge. 

PACS numbers: 04.20.Jb 
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In Einstein-Maxwell (EM) theory the unique family of stationary asymptotically flat black holes comprises the 
rotating Kerr-Newman (KN) and Kerr black holes and the static Reissner-Nordstr0m (RN) and Schwarzschild black 
holes. EM black holes are uniquely determined by their mass M, their angular momentum J, their electric charge Q, 
5^ | and their magnetic charge P, i.e. EM black holes have "no hair" 0,0- 

&)[)■ The EM "no-hair" theorem does not readily generalize to theories with non-Abelian gauge fields coupled to gravity 
0,0. The hairy black hole solutions of SU(2) Einstein- Yang- Mills (EYM) theory possess non-trivial magnetic fields 
outside their regular event horizon, but carry no magnetic charge 0, 0,0, 00- Their magnetic fields are characterized 
by the winding number n and by the node number k of the gauge field functions. In the static limit, the solutions 
with winding number n — 1 are spherically symmetric, whereas the solutions with winding number n > 1 possess only 
axial symmetry, showing that Israel's theorem does not generalize to non-Abelian theories, either 

In many unified theories, including string theory, dilatons appear. When a dilaton is coupled to EM theory, this has 
profound consequences for the black hole solutions. Although uncharged Einstein-Maxwell-dilaton (EMD) black holes 
correspond to the EM black holes, since the source term for the dilaton vanishes, charged EMD black hole solutions 
possess a non-trivial dilaton field, giving rise to an additional charge, the dilaton charge D, and to qualitatively new 
features of the black holes. Charged static EMD black hole solutions, for instance, exist for arbitrarily small horizon 
size 0, and the surface gravity of 'extremal' solutions depends in an essential way on the dilaton coupling constant 7. 
Rotating EMD black holes, known exactly only for Kaluza-Klein (KK) theory with 7 = \/3 no longer possess 

the gyromagnetic ratio g = 2 the value of KN black holes. Extremal charged rotating EMD black holes can 
possess non-zero angular momentum, while their event horizon has zero angular velocity |10| . 

Here we consider rotating black hole solutions of SU(2) Einstein- Yang-Mills-dilaton (EYMD) theory 
For fixed dilaton coupling constant 7 and winding number n, the black hole solutions form sequences, which with 
increasing node number k, tend to limiting solutions. The black hole solutions of a given sequence carry no non- 
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Abelian magnetic charge, but generically they carry a small non-Abelian electric charge [|| 0, El- The limiting 
solutions, in contrast, correspond to embedded Abelian black hole solutions, which carry no electric charge, but a 
magnetic charge, equal to the winding number n. 

We investigate the physical properties of these black hole solutions. In particular, we consider their global charges, 
obtained from an asymptotic expansion of the fields. The asymptotic expansion of the gauge field functions generically 
involves non-integer powers of the radial coordinate, depending on the non-Abelian electric charge of the black hole 
solutions 0. In the presence of the dilaton, there is, however, also a special set of solutions, whose non-Abelian 
electric charge vanishes Q- For these, the asymptotic expansion involves only integer powers of the radial coordinate. 
We further consider the horizon properties of the non-Abelian black hole solutions, such as their horizon area A, their 
surface gravity K sg , their horizon curvature K and horizon topology, and we introduce a topological horizon charge 

v 0.GJ. 

For the rotating non-Abelian black holes the zeroth law of black hole mechanics holds [a, Q , as we ^ as a generalized 
first law ^{|. The non-Abelian black holes further satisfy the mass formula 0, 

M = 2TS + 20J + — + 2ip cl Q , (1) 

7 

where T denotes the temperature of the black holes, S their entropy, f2 their horizon angular velocity, and ij) e i 
their horizon electrostatic potential. This non-Abelian mass formula is similar to the Smarr formula 16]. However, 
instead of the magnetic charge P, the dilaton charge D enters. This is crucial, since the genuinely non-Abelian black 
hole solutions carry magnetic fields, but no magnetic charge. Thus the dilaton charge term takes into account the 
contribution to the total mass from the magnetic fields outside the horizon. This mass formula holds for all non- 
perturbatively known black hole solutions of SU(2) EYMD theory p|, 0, Il2l . including the rotating generalizations 
of the static non-spherically symmetric non-Abelian black hole solutions p£. It also holds for embedded Abelian 
solutions HE3- 

Concerning the uniqueness conjecture for EYMD black holes, it is not sufficient to simply replace the magnetic 
charge P by the dilaton charge D. Black holes in SU(2) EYMD theory are not uniquely characterized by their mass 
M, their angular momentum J, their dilaton charge D, and their non-Abelian electric charge Q. Adding as an 
additional charge a topological charge N, however, a new uniqueness conjecture can be formulated p|: Black holes in 
SU(2) EYMD theory are uniquely determined by their mass M, their angular momentum J, their dilaton charge D, 
their non-Abelian electric charge Q, and their topological charge N. 

In section II we recall the SU(2) EYMD action and the equations of motion. We discuss the stationary ansatz for 
the metric and the gauge and dilaton fields, and we present the boundary conditions. In section III we address the 
properties of the black hole solutions. We briefly present the asymptotic expansion at infinity and at the horizon, 
from which we obtain the global charges and the horizon properties, as well as the proof of the mass formula. Our 
numerical results are discussed in section IV. In section V we present our conclusions. Appendices A and B give 
details of the expansion at infinity and at the horizon, respectively. In Appendix C we discuss the energy conditions 
and algebraic type of the stress-energy tensor in SU(2) EYMD theory. 



II. NON-ABELIAN BLACK HOLES 



After recalling the SU(2) EYMD action and the general set of equations of motion, we discuss the ansatz for 
the stationary non-Abelian black hole solutions, employed for the metric, the dilaton field and the gauge field 0. 
The ansatz for the metric represents the stationary axially symmetric Lewis-Papapetrou metric |l8j in isotropic 
coordinates. The ansatz for the gau ge field includes an arbitrary winding number n and satisfies the Ricci 

circularity and Frobenius conditions 18]. As implied by the boundary conditions, the stationary axially symmetric 
black hole solutions are asymptotically flat, and possess a regular event horizon. 

A. SU(2) EYMD Action 

We consider the SU(2) Einstein- Yang-Mills-dilaton (EYMD) action 

5 = + (2) 

with scalar curvature R and matter Lagrangian Lm given by 

Lm = —B^V* - ie 2K *Tr(^FH , (3) 
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with dilaton field <&, field strength tensor F M „ = d^A^ — d v A^ + ie [A^, A v \, gauge field A^ = A°r a /2, and Newton's 
constant G, dilaton coupling constant k, and Yang-Mills coupling constant e. 

Variation of the action with respect to the metric and the matter fields leads, respectively, to the Einstein equations 

G^v — R^v — -g^R — 8ttGT^ (4) 

with stress-energy tensor 

.dL M 



d^dv® - \g^d a $d a <S> + 2e 2K *Tr(F MQ F„ /35 ^ - lg^F a0 F al3 ) , (5) 



and the matter field equations, 

1 



(V^gd^) = K e 2K *Tr (F fiU F flv ) , (6) 



--D^^e^F^) = , 



(7) 



where A, = cL + ie [A,, •]. 



B. Stationary Axially Symmetric Ansatz 

The system of partial differential equations, (QJ, 10, and (7J|, is highly non- linear and complicated. In order to 
generate solutions to these equations, one profits from the use of symmetries, simplifying the equations. 

Here we consider black hole solutions, which are both stationary and axially symmetric. We therefore impose on 
the spacetime the presence of two commuting Killing vector fields, £ (asymptotically timelike) and r\ (asymptotically 
spacelike). Since the Killing vector fields commute, we may adopt a system of adapted coordinates, say {t,r, 0, ip}, 
such that 

t = d t , V = d v . (8) 

In these coordinates the metric is independent of t and if. We also assume that the symmetry axis of the spacetime, 
the set of points where r\ = 0, is regular, and satisfies the elementary flatness condition 

^#^ = 1> X = ^ • (9) 

Apart from the symmetry requirement on the metric [C^g = C^g = 0, i.e., g^ v = g^fr^O)), we impose that the 
matter fields are also symmetric under the spacetime transformations generated by £ and r/. 
This implies for the dilaton field 

£ c $ = = , (10) 

so $ depends on r and 9 only. 

For the gauge potential A = A^dx^, the concept of generalised symmetry 0,|2(j requires 

(C^A), = D,Wt , 

(C V A)^ = D^W V , (11) 
where W% and W v are two compensating su(2)-valued functions satisfying 

C e W v - C v W s + ie[W^Wr,] =0 . (12) 
Performing a gauge transformation to set W% — 0, leaves A and W v independent of t. 
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To further simplify the system of equations, one can impose that the Killing fields generate orthogonal 2-surfaces 
(Frobenius condition) 



£AriAd€ = £AT)AdT) = 0, 



(13) 



where £ and r\ are considered as 1-forms. By virtue of the circularity theorem [2l| such a condition may be written 
in terms of the Ricci tensor 



(14) 



£ A T7 A R(£) = £ A rj A R(r)) = , 
where (R(v))p = R^v" '. 

The metric can then be written in the Lewis-Papapetrou form |22j | , which in isotropic coordinates reads 

ds 2 = -fdt 2 + - \dr 2 + r 2 d6 2 } + sin 2 Or 2 - \dip - -dt] * , (15) 
/ / L r J 

where /, m, I and u> are functions of r and 9 only. Note the change of sign of the function ui with respect to Ref. [||. 
The z-axis represents the symmetry axis. The regularity condition along the z-axis © requires 

m\e=o,7T = i|e=o,7r • (16) 

The event horizon of stationary black hole solutions resides at a surface of constant radial coordinate, r = rn, and 
is characterized by the condition /(rn, 0) = 6]. The Killing vector field 



x = £ + — 



(17) 



is orthogonal to and null on the horizon [23j. The ergosphere, defined as the region in which is positive, is 

bounded by the event horizon and by the surface where 



/ + sin 2 6 j lj 2 = . 



(18) 



Due to the Einstein equations the circularity conditions (|14|) have consequences on the matter content, namely, 

CAT 7 AT(e) = CAT 7 AT(r 7 ) = 0, (19) 

with (T(w)) A1 = T^ v v v . However, contrary to the case of Abelian fields, the conditions (|19fl are not just a consequence 
of the symmetry requirements, but they give rise to additional restrictions on the form of the gauge potential. 

For the gauge fields we employ a generalized ansatz yj, which trivially fulfils both the symmetry constraints i|ll|) 
and l|12f> and the circularity conditions 119f) . 



= mt + A v {dtp - -dt) + ( ^-dr + (1 - H 2 )d9 ) 



(20) 



2e 2e 



(21) 



A, 



-n sin 6 



(22) 



Here the symbols r™, Tg and r™ denote the dot products of the Cartesian vector of Pauli matrices, f = (t x ,t v ,t z ), 
with the spatial unit vectors 



e r = (sm a cos nip, smOsm nip, cos a) 
eg 1 = (cos 8 cos nip, cos 6 sin nip, — sin t 
e" = (— smnip, cos nip, 0) , 



(23) 



respectively. Since the gauge fields wind n times around, while the azimuthal angle ip covers the full trigonometric 
circle once, we refer to the integer n as the winding number of the solutions. The gauge field functions Bi and Hi 
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depend on the coordinates r and 9 only. (For n — 1 and k = 0, the previously employed stationary ansatz Q is 
obtained, whereas for uj = B\ = B 2 — the static axially symmetric ansatz is recovered including the static 
spherically symmetric ansatz [lj], for n — 1 and Hi = H3 = 0, H2 = -H4 = w(r) and $ = 
The ansatz is form-invariant under Abelian gauge transformations U 



[/ = exp ( VTM) I . (24) 



With respect to this residual gauge degree of freedom we choose the gauge fixing condition rd r H\ — d$Hi = 0-0] 
For the gauge field ansatz, Eqs. l]2U p -l]22 [l . the compensating matrix W v is given by 



W v = n^- . (25) 



We note, that by employing the gauge transformation U |2J 



U = exp ( -T z mp ) , (26) 



one can choose W' n = |19|, |25| , leading to the gauge field A'^ 



1"' 



A'dx* = 9'dt + A 1 (dip - -dt) + ( —dr + (1 - H 2 )d6 ) ^, (27) 
F ^ r \ r J 2e 



= Si (sin (9^ + 008(9-^") + 5 2 fcos 9^- - sin 0^) - n^-- , (28) 
V 2e 2e/ V 2e 2e/ 2e r y ' 



A — —nsmt 



H 3 ( sin 0^- + cos0^-) + (1 - HA fcos 0^ - sin 9^-)} - n^-. (29) 
V 2e 2e/ V 2e 2e/J 2e 



We further note, that when imposing restricted circularity conditions by requiring F(£,7j) — *F(£,r)) — 0, only 
Abelian solutions are possible if asymptotic flatness is assumed pfij . 

C. Boundary Conditions 

Dimensionless Quantities 

For notational simplicity we introduce the dimensionless coordinate x, 

x = -^=r , (30) 



the dimensionless electric gauge field functions B\ and £? 2 , 



ft.^ft, 5 2 = ^ B2 , ,31) 
e e 



the dimensionless dilaton function 



<t> = V4^G$ , (32) 

and the dimensionless dilaton coupling constant 7, 

For 7 = 1 contact with the low energy effective action of string theory is made, whereas in the limit 7 — > the 
dilaton decouples and EYM theory is obtained. 

Boundary conditions at infinity 
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To obtain asymptotically flat solutions, we impose on the metric functions the boundary conditions at infinity 

f\x— oo — m\x= oc — ^|a;— oo — 1 ; oo — ■ (34) 

For the dilaton function we choose 

<£|*=oo = , (35) 

since any finite value of the dilaton field at infinity can always be transformed to zero via 4> — > </>— 4>(oo), x — > xe -7 ^ 00 '. 
We further impose, that the two electric gauge field functions Bi vanish asymptotically, 

-Bl|x=oo = B2\x=ao = • 

For magnetically neutral solutions, the gauge field functions Hi have to satisfy 

H\\ x=00 — Hz\ x=00 = , H 2 \ x =oo = H4\ x=OQ = (— l) fe , (36) 

where the node number k is defined as the number of nodes of the functions H2 and H4 along the positive (or negative) 
z-axis Note, that for each node number there is a degenerate solution with i?2U=oo = #4 1 3=00 = — ( — l) fc ) 

related by the large gauge transformation U = ir™. Under this transformation the functions transform according to 

Hi — > —Hi , H2 — ► — H 2 , H 3 — > +H 3 , H4 — > —Hi , 

B 1 -*+B 1 , B 2 -> -B 2 -2nsm6- . 

x 

Boundary conditions at the horizon 

The event horizon of stationary black hole solutions resides at a surface of constant radial coordinate, x = xn, and 
is characterized by the condition f(xn, 0) = |6j. 

Regularity at the horizon then requires the following boundary conditions for the metric functions 

f\ x =x H = m\ x=XH = l\ x=XH = , oj\x=x a = = const. , (37) 

for the dilaton function 

»«^U=xh = . (38) 

and for the 'magnetic' gauge field functions 

Hi\x=x H — j d x H 2 \ x=XH = , d x H 3 \ x=XH = , d x H4 : \ x=XH = , (39) 

with the gauge condition dgHi = taken into account 

As discussed in 6], the boundary conditions for the 'electric' part of the gauge potential are obtained from the 
requirement that the electro-static potential "J = X^-^-n is constant at the horizon [27| . 

§(r H )=X ,1 ^| r=rH =* H ■ (40) 

Defining the dimensionless electrostatic potential ip, 



rf> = Vi^G* , (41) 

and the dimensionless horizon angular velocity Q, 

n = ^, (42) 

xn 

this yields the boundary conditions 

Bi\x=xn = nQcos9 , B 2 \ x=Xll = -nfisinfl . (43) 
With these boundary conditions, the horizon electrostatic potential reads 

^ H = nO^ = Vciy , (44) 
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defining ip e \ for the non-Abelian mass formula. 

Boundary conditions along the axes 

The boundary conditions along the p- and z-axis (8 = ir/2 and 9 = 0) are determined by the symmetries. They are 
given by 

def\e=o = dem\ e =o = d e l\e=o = d g uj\e=o = , 

<9e/|e=f = d s m\ B= iL = dgl\ s= % = d e uj\ e =% = , (45) 
de<p\ e= a = , 

<9 e 0| 9=f = , (46) 

-82(0=0 = , dgBi\g= = , 
-ffl|e=0 = H 3 \g= = , dgH 2 \g=o = dgH4\g=o = , 

Bi\ 0=i =0, dgB 2 \ e=i =0, 

Hi\e=i = Ha\e=i - , dgH 2 \ e=i = d e H 4 \ e ^ = . (47) 

In addition, regularity on the z-axis requires condition (|16J) for the metric functions to be satisfied, and regularity of 
the energy density on the z-axis requires 

H 2 \g=o - ff 4 |e=o ■ (48) 

III. BLACK HOLE PROPERTIES 

We derive the properties of the stationary axially symmetric black holes from the expansions of their gauge field and 
matter functions at infinity and at the horizon. The expansion at infinity yields the global charges of the black holes 
and their magnetic moments. Generically the gauge field functions of the black hole solutions show a non-integer power 
fall-off asymptotically, with the exponents determined by the non-Abelian electric charge Q. The expansion at the 
horizon yields the horizon properties, such as the area parameter, the surface gravity, and the horizon deformation. 
We also introduce a topological charge of the horizon. We then give a detailed account of the non-Abelian mass 
formula. 

A. Global Charges 

The mass M. and the angular momentum J of the black hole solutions are obtained from the metric components 
gtt and gt v , respectively. The asymptotic expansion for the metric function /, 

x \x z 
yields for the dimensionless mass M, 

M = ~ lim x 2 d x f , (49) 

Z X — >QC 

where 



M = ^pM . (50) 
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Likewise, the asymptotic expansions for the metric function u, 

2J „ ( 1 

X z \X' 

yields for the dimcnsionless angular momentum J, 



where 



J = — lim x 2 uj , (51) 

2 x — >oc 



Air 

J = ^J. (52) 



The asymptotic expansion for the dilaton function 4>, 

D iQ\ (i 



x 2x 2 

yields the dimensionless dilaton charge D, 

D = lim x 2 d x (j) , (53) 

x — >oo 

related to the dilaton charge T> via 

V = -D . (54) 
e 

The asymptotic expansion of the gauge field yields the global non-Abelian electromagnetic charges, Q and V . The 
asymptotic expansion of the electric gauge field functions B\ and B 2 , 

B^Qmal+oa), B 2 = - { -it^ + (\ 

x \ x * J x \ r 



yield the dimensionless non-Abelian electric charge Q, 

Q = - lim x 2 d x {cosOBi - (-l) fe sm9B 2 ) , (55) 

x— ►OO 

where 

Q = — . (56) 
e 

The boundary conditions of the magnetic gauge field functions guarantee that the dimensionless non-Abelian magnetic 
charge P vanishes. 

The non-Abelian global charges Q and P appear to be gauge dependent. In particular, the definition of Q corre- 
sponds to rotating to a gauge, where the gauge field component ip = V '^G^ asymptotically only has a r z component 


y x 2 

Identifying the global non-Abelian electric charge Q in this way, corresponds to the usual choice 0, 12H l2Sj . 

We note, that the modulus of the non-Abelian electric charge, \Q\, corresponds to a gauge invariant definition for 
the non-Abelian electric charge given in Ref. |2*9| . 



where *F represents the dual field strength tensor, and the integral is evaluated at spatial infinity. 
Likewise, a gauge invariant definition of the non-Abelian magnetic charge is given in Ref. [2j| , 



where again the integral is evaluated at spatial infinity, yielding P = 0. 

Note that the lowest order terms in the expansions, as needed for the expressions of the global charges M, J, D, Q 
and P, do not involve non- integer powers. 
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B. Non-integer power fall-off 

Here we present the lowest order terms of the expansion of the 'magnetic' gauge field functions H± - H4 for winding 
number n = 1 and n — 3. (For details see Appendix A). Since for n = 2 the lowest order expansion contains non- 
analytic terms like \og(x)/x 2 already in the static limit |30j |. we refrained from an analysis of the rotating solutions in 
this case. 

The asymptotic expansion yields for n — 1 



Hi 

Hi 
H 3 



2C5 8C4 -j(0-i) _ 2C 2 C 3 (a + 3) i fa+1) 
x 2 (3-1 (a + 5)Q 2 



Sill (-) CI IS (-) + O I -^r 



(Q- - (-l) fe C 3 ar5(«-i)J sin0cos0 + oQ 
Hi = {-if + C^x-^-^ + ({-l) k — -CaX-^-V ) sin 2 6 + [ -) , (59) 







) sin 2 e 









where C, are dimensionless constants. 

Here a and (3 determine the non-integer fall-off of the gauge field functions, 



a = V9 - 4Q 2 , ft = a/25 - 4Q 2 . (60) 
For n = 3 the lowest order terms of the expansion involve only integer powers, 

H x = 5hin2<9 + o(4j 

% A \x z 

H 2 = (-l) k + ^.cos26 + o(J^j , 

H 3 = — sm cos H — ; sin cos 6/ + o — 

a; 2x z \x z 

h 4 = (-D* + (-1)^ sin2 fl + % + (_!)* + ^ - - ^ sin2 e 

x x z 2x z 

(61) 

but higher order terms contain non-integer powers involving e = \/49 — 4Q 2 . 

To determine, whether the generic non-integer fall-off is a physical property of the solutions or a gauge artefact, 
we consider the asymptotic behaviour of the gauge invariant quantities Tr (F^F^ V ) and Tr {F^„* F^ v ). Inserting the 
expansions Eqs. (|59|l and 161|) we find 

KiF» v Fn = -% + m -l D)Q \ o(\) , „ = landn = 3, (62) 

X X \ X J 

and 

Tr (F^Fn = l^ C os0 + o ^) , „=1, 

Tr^/F-) = i^cos0 + O (-l) , n = 3. 

Although lowest order terms do not contain the non-integer powers, they occur in the higher order terms, indicating 
that the non-integer power decay cannot be removed by a gauge transformation. 
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C. Magnetic Moment 

In Abclian gauge theory, the dimensionless magnetic field Ba of a magnetic dipole with dipole moment ft A = fife z 
is given by 

A 

B A = f\ (2cos6»e r + sinfle e ) . 

To investigate the magnetic moment for the non- Abclian black holes we work in the gauge where asymptotically 
Aq = + o(l jx). From the asymptotic expansion of the gauge field functions, given in Appendix A, we find for 
the magnetic field for winding number n = 1 solutions 

B f = -^§2cos0^+&^ , 

„ C 2 T z T 

= H T \ , (63) 

where the functions bf,bg,b0 involve non- integer powers of x, and contain the leading terms in the expansion for 
Q 7^ 0. Since the leading terms decay slower than a; -3 for Q ^ we cannot extract the magnetic moment in general. 

In the following we give a suggestion for the definition of the magnetic moment. The asymptotic form of Ao suggests 
to interpret ?f as 'electric charge operator' Q e su(2) in the given gauge. We then observe that the projection of Bp 
in Q-dircction, 



B^> = -^2cos9 + o(l/x a ) , 
Bf ] = -^-sm6 + o(l/x 3 ) , 



Bf = o(l/x 3 ) (64) 
corresponds asymptotically to the magnetic field of a magnetic dipole ft = —C 2 e z with magnitude 

iMl = \Ci\ ■ (65) 
We note that ft is invariant under time independent gauge transformations, since such a gauge transformation rotates 
the 'electric charge operator' Q and the magnetic field Bp in the same way and leaves the projection B^ invariant. 

A special case arises for n = 1 and Q = 0, when the functions bf,bg,b^ are of order 0(x~ 3 ). In this case the 
asymptotic form of the magnetic field reads 

B f = -^2cos0^ + ^2sin0^ + O(l/a; 4 ) , 
x 6 2 x 6 2 



B = -9l T A + o(l/x*), (66) 



x° 



Comparison with the magnetic field of a magnetic dipole ft A in Abclian gauge theory, 



x° 

= — (yU^(2cos6>e r + sin flee) + /U^(2sin^e r — cos flee) + ^e v ) 



X' 

suggests the identification of the su(2) valued magnetic moment 



= -C 2 y . (67) 
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In the static case C3 = — C 2 , reflecting the spherical symmetry of the solution. For the stationary rotating solutions 
the spherical symmetry is broken to axial symmetry, corresponding to C3 =/= — C 2 . 

For winding number n = 3 the asymptotic expansion of the gauge field functions yields for the magnetic field 



leading again to /2 — —C 2 e z - 



B A 



.9l2cos6^ + 0(l/x A ) , 
x i 2 

C 2 



smf 



0(l/x 4 



= 0(l/x 4 ) 



D. Expansion at the Horizon 



(68) 



Expanding the metric and gauge field functions at the horizon in powers of 



yields to lowest order 



me 

m{5,d 
1(5,9 
uj(S,9 
<t>{8,0 

B 2 (S,9 

H 2 {8,9 
H 3 (S,9 
H 4 (S,9 



S= — -l 

XH 



S 2 f 2 (l-S) + 0(S 4 ) , 
6 2 m 2 {l-35) + 0{5 4 ) , 
5 2 l 2 (l-35) + 0(S 4 ) , 
wh(1 + 6) + 0(S 2 ) , 
0o + 0(S 2 ) , 

n— cos9 + 0{5 2 ) , 
-n—sm0 + O(6 2 ) , 

H 20 + 0(5 2 ) , 
#30 + 0(6 2 ) , 
ff 40 + 0{8 2 ) , 



(69) 



(70) 



The expansion coefficients f 2 , m 2l l 2 , 4>o, Hn, H 2 o, H30, H^q are functions of the variable 9. Among these coefficients 
the following relations hold, 



_ d g m 2 _ 2 def 2 
m 2 " f 2 



H\\ — dgH 20 . 

Further details of the expansion at the horizon are given in Appendix B. 



(71) 
(72) 



E. Horizon Properties 

Let us now discuss the horizon properties of the SU(2) EYMD black hole solutions. The first quantity of interest 
is the area of the black hole horizon. The dimensionless area A, given by 



A = 2ir d9 sin 



>\lhjn 2 



(73) 



defines the area parameter xa via |2£ 



A = 4ttx a , 



(74) 
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and the dimensionless entropy S of the black hole 



A 
1 



The surface gravity of the black hole solutions is obtained from [2S 

1 



-(v^)(vy) , 



(75) 



(76) 



with Killing vector x, Eq. fTTjl . Inserting the expansion at the horizon, Eqs. J2U, yields the dimensionless surface 
gravity 



(77) 



As seen from Eq. (|71(l . K sg is indeed constant on the horizon, as required by the zeroth law of black hole mechanics. 
The dimensionless temperature T of the black hole is proportional to the surface gravity, 



2tt 



(78) 



To obtain a measure for the deformation of the horizon we compare the dimensionless circumference of the horizon 
along the equator, L c , with the dimensionless circumference of the horizon along a great circle passing through the 
poles, L pi 




L p = 2 



x—xn ,9= 7r/2 



dO 4 / — x 



x— xh ,ip— const. 



We obtain further information about the horizon deformation by considering its Gaussian curvature K, 



K{9) = 



.92 



, 92 = geeg vv 



We determine the topology of the horizon, by computing its Euler characteristic xe. 



XE 



1 

2^ 



Using the expansion Eqs. (|70|l in the integrand, 



K\f92\x=xn — ~~qq 



cos9\l — 

m,2 



;2 

'2 



sin0 d 



2y/h d9 V m 2 



(79) 



(80) 



(81) 



(82) 



and the fact that li = mi on the z-axis, this yields for the Euler number \e = 2, indicating that the horizon has the 
topology of a 2-sphere. 

To find the horizon electric charge Q\ and the horizon magnetic charge Pa, on e can evaluate the integrals, Eqs. (|57|> 



and l|58|) . representing dimensionless gauge invariant quantities |29|. Whereas such a definition of the non-Abelian 
horizon charges appears adequate in the static case, it appears problematic in the stationary case. For embedded 
Abelian solutions, for instance, one obtains a horizon electric charge, which differs from the global electric charge. The 
reason is that when evaluating the horizon electric charge according to this prescription, one is taking the absolute 
value of the dual field strength tensor. Thus one does not allow for the cancellation, present in a purely Abelian 
theory, when the integral involves the dual field strength tensor. 

We finally introduce a topological horizon charge TV, suggested by Ashtekar [l4|. For that purpose we consider the 
su(2)-algebra valued 2-form Fh, corresponding to the pull-back of the Yang-Mills field strength to the horizon H , 



Fh = F ev \nd9 f\dip 



(83) 



Its dual *Fr on the horizon 



*F H = g 2 1/2 F g , 



(84) 
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represents an su(2)-algebra valued function on the horizon. We decompose *Fh as 

*F H = *F£ (cosn^ + sinn^) + , (85) 

with norm | *-Fh|, 



I^h| = V(^h) +(*F£) . (86) 
Next we define a normalised su(2) valued function on the horizon by 

= sin 9 I cosnip— + s'mmp— ) + cosG— , (87) 



|*F H | V 2 T 2l 2 

where 

smQ = *F^/\*F B \ , coaQ = *F£/\*F s \ . 
a represents a map from the 2-sphere of the horizon to the 2-sphere in the Lie-algebra su(2). The degree of the map, 



N = ^- [ \e abc a a da b A da c , 



4^ 7 H 2 



defines the topological horizon charge N, 



If n 

N= — (sin ®nd e ®)d6d<p = —- (cos 9\g=„ - cos Q\e=o) • (89) 
47T J H 2 

To evaluate the topological horizon charge N for the stationary non-Abelian black hole solutions, we express cos Q 
and sin O in terms of the functions Hi . This yields 

. ^ sin 8G r + cos OGg „ cos9G r — sm9Gg 

suit) = ; , cost) = ; , (90) 

where 

G r = - [H 3 . e + H 3 cot6 + H 2 H 4 - 1] - , 

e 

G e = [i/4.0 + cot0(if 4 - H 2 ) - H 2 H 3 ] - . 

e 

Symmetry and boundary conditions imply G r \g =T! — G r \g = o and Gg\g =T! — Gg\g=o = 0, therefore 



N = 



G r 



nsign(G r | 9= o) , (91) 



0=0 



relating the winding number n to the topological horizon charge N. We note that solutions with winding number — n 
are gauge equivalent to solutions with winding number n, but carry opposite electric charge |3lj . 

F. Non-Abelian Mass Formula 
The non-Abelian mass formula, Eq. 0), 

M = 2TS + 2n,J + — + 2^ H Q 

1 

is a generalization of the non-Abelian mass formula M = 2TS + -D/7, obtained previously for static axially symmetric 
non-Abelian solutions 0, |2|J, which generically carry no non-Abelian charges. Rotating non-Abelian black hole 
solutions in general do carry non-Abelian electric charge, but they do not carry non-Abelian magnetic charge. Since 
they nevertheless carry non-trivial non-Abelian magnetic fields, the Smarr formula |l6|. 

M = 2TS + 2dJ + ^cl<2 + ^magi 3 , (92) 
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(where ip mag represents the horizon magnetic potential) is bound to fail: The contribution to the mass from the non- 
Abelian magnetic fields outside the horizon would not be included. In constrast, in the non-Abelian mass formula 

i ffl l. this magnetic field contribution to the mass is contained in the dilaton term — . 

7 

To derive the non-Abelian mass formula, let us begin by considering the general definition of the local mass A^iocai 
and local angular momentum j7i OC ai, (see e. g. |23j ) 



M 



local 



— I- 

8ttGJ s 



(93) 



167rG J S 2 2 



(94) 



for a given 2-sphcrc S 2 . By taking the 2-sphcre to infinity, and employing Stokes theorem, the global charges M. and 
J are obtained j2^|, 



M = 



1 



4ttG 



1 



8ttG 7h 2 



1 

—e 



pv pa 



\l p i a dx ii dx v 



(95) 



J = ~7^ f R^rfdV+-±- I \ep, vpa Wrfdx»dx v 
8ttG J s ^ 16ttG J h 2 



(96) 



Here £ denotes an asymptotically flat hypersurface bounded by the horizon H, = (1, 0, 0, —uj/r)/\/J is normal to 
E with ripn^ = —1, and dV is the natural volume element on £ (dV = m^l/f 3 r 2 sm8drd8d<p), and consequently, 



*b i R ^ dV =-^G j^^drd0d V 



(97) 



Defining the horizon mass of the black hole Mn = -Miocaiifn) and its horizon angular momentum J7k = j7iocai(?"H) 
one obtains with \ — £ + (^h/'"h) V the horizon mass formula 



M h = 2TS + 2^Jh , 

where T and S denote temperature and entropy. The global mass is then given by the relation 



l — J^R° V=gdrdOdcp . 



M = 2TS + 2— Jh - . ,., 

r H 47rG j s 

Let us next replace the integrand in (|99|l by first making use of the Einstein equations, 



(98) 



(99) 



R% = 16vrGe 2K *Tr ( F 0fJ _F 0fJ - - -F^F"" 



(100) 



and then of the dilaton field equation, 



i?° n = 8vrG 



1 1 



We now evaluate the integral involving the dilaton d'Alembertian 
1 



gd^) + 2e 2K<I 'Tr (F pF 0p ) 



-.dp (y/=g&*$) 



gdrdOdip = [ [d r (^d r ^>) + d g (y/=g~d e $)] drdOdy , 
is 



(101) 



(102) 



with the divergence theorem. Making use of the boundary conditions at the horizon (9 r <f> = 0) and at infinity 
(lmv^oo \/—gd r Q = 2? sin 6*), yields for the first term in l|102|) 



2tt / V^(jd r <S>\ r r z™d9 = AttT> , 
'o 



(103) 
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with the dilaton charge T>. The second term vanishes, since y/—g vanishes at 8 = and 6 = %. The relation for the 
global mass then reads 



M = 2TS + 2— Jh + — V - 4 / e 2K *Tr[F „F°' i ]</=g~drddd<p 



(104) 



Making again use of the Einstein equations, we next replace the horizon angular momentum J7k by the global 
angular momentum J pflj . 



J = Jn + 2 / e 2K *Tr(F w F »)^jdrd8dcp , 



and obtain 



M - 2TS - 2— J - —V = X : 



2 2K *Tr 



X, 



0/i 



F w 1 X ^ 1 



-gdrdOdtp , 



defining the integral X. 

To evaluate X, Eq. (|106|> . we make use of the symmetry relations, Eqs. (|ll|l 20], 



i^o = D^A , F w = X> M (A v - W,,) 



The integral then reads 



1 = 4 



e 2K *Tr 



\d, (a + ^ - ^)) } 



' '—gdrdOdtp 



(105) 



(106) 



(107) 



(108) 



Adding zero to the above integral, in the form of the gauge field equation of motion for the zero component, the 
covariant derivative then acts on all functions in the integrand, 



I = 4 jf Tr | (A + ^ (4, - W„) ) e 2K *^ V^ff} 



drdfldp . (109) 
Since the trace of a commutator vanishes, we replace the gauge covariant derivative by the partial derivative, 

drdOdip , (110) 



X = 4 I Tr 



^ \ ( A + ^ {Av - W v ) ) e 2 ^F^V~g 



and employ the divergence theorem. The 0-term vanishes, since \J—g vanishes at 6 = and 9 = n, and the (/3-term 
vanishes, since the integrands at (p = and </? = 2-7T coincide, thus we are left with 



X = 4 / Tr 



A + — (A v -W v )]e 



rn 



2k* rjiOr 



dOdtp 



(111) 



We next make use of the fact that the electrostatic potential \P|h is constant at the horizon (see Eqs. (|40|l and 144|) '). 



= [A + ^A V 



r H 2 



(112) 



Hence the integrand vanishes at the horizon, and the only contribution to X comes from infinity. 
At infinity the asymptotic expansion yields to lowest order 



F 0r V^g = --sine I cos0^- - (-l) ft sin 



Q 



A = o(l) , 

A v = --(1 -(-!)*) Bin^ 



Hence Aq does not contribute to the integral, and A v contributes for odd node number k. The integral X, 



iwnQ wh 



'h 



(113) 



(114) 
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is therefore independent of k, and the mass formula becomes 

Mtt AlT 

M-2TS-2—J P = 87r*oiQ. (115) 

r H k 

To complete the proof of the mass formula (QJ for genuinely non-Abelian black holes, subject to the above ansatz 
and boundary conditions, we now return to dimensionless variables. Noting that 



TS = ^-TS , (116) 

- = -vh n > 

and recalling the definition of ipe\, Eq. 1441) . the non-Abelian mass formula is obtained, 

M = 2TS + 2flJ + — + 2ifj cl Q . 

7 

The mass formula also holds for embedded Abelian black holes 0, • 

IV. NUMERICAL RESULTS 

We solve the set of eleven coupled non-linear elliptic partial differential equations numerically [12, subject to 
the above boundary conditions, employing compactificd dimensionless coordinates, x = 1 — (x-r/x). The numerical 
calculations, based on the Newton- Raphson method, are performed with help of the program FIDISOL 32] . The 
equations are discretized on a non-equidistant grid in x and 9. Typical grids used have sizes 100 x 20, covering the 
integration region < x < 1 and < 9 < n/2. (See 0,0, H(| and [33 for further details on the numerical procedure.) 

For a given dilaton coupling constant 7, the rotating non-Abelian black hole solutions then depend on the continuous 
parameters xn and uh, representing the isotropic horizon radius and the metric function to at the horizon, respectively, 
while their ratio whAh corresponds to the rotational velocity of the horizon. The solutions further depend on the 
discrete parameters n and fc, representing the winding number and the node number, respectively. 

To construct a rotating non-Abelian black hole solution of SU(2) EYMD theory with dilaton coupling constant 7 
and parameters ih, wh, n>, and k, we either start from the static SU(2) EYMD black hole solution possessing the 
same set of parameters except for wh = 0, and then increase wh or we start from the corresponding rotating SU(2) 
EYM black hole solution and then increase 7 |g. 



A. n = 1 Black Holes 



Rotating non-Abelian black hole solutions with winding number n = 1 have been studied before in EYM theory 
. Here we study the influence of the presence of the dilaton on the properties of these black hole solutions, with 
particular emphasis on the Q = black hole solutions. 

Global charges 

Many features of the n = 1 EYMD black hole solutions agree with those of the n — 1 EYM solutions, studied before 
. In particular, when increasing wh from zero, while keeping xh fixed, a lower branch of black hole solutions forms 
and extends up to a maximal value cj^ ax (xH,7), where an upper branch bends backwards towards = 0. This is 
seen in Figs, la-d, where the mass M, the angular momentum per unit mass a = J/M, the relative dilaton charge 
-D/7, and the non-Abelian electric charge Q are shown as functions of the parameter ton for black holes with horizon 
radius ih = 0.1, winding number n = 1, and node number k = 1 for five values of the dilaton coupling constant, 
7 = 0,0.5, 1, a/3 and 3. For 7 = the relative dilaton charge -D/7 is extracted from the mass formula Eq. 

The mass M (Fig. la) and the angular momentum per unit mass a (Fig. lb) of the non-Abelian solutions increase 
monotonically along both branches, and diverge with ui^ 1 in the limit - > along the upper branch. Both mass 
and angular momentum become (almost) independent of the dilaton coupling constant along the upper branch. The 
relative dilaton charge -D/7 (Fig. lc) decreases monotonically along both branches, approaching zero in the limit 
lu-r — * 0. (The curves are discontinued along the upper branch when the numerical procedure no longer provides 
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high accuracy.) Also the dilaton charge approaches its limiting value along the upper branch with a slope, (almost) 
independent of the dilaton coupling constant. 

The non-Abelian electric charge Q (Fig. Id) increases monotonically along both branches, when 7 < 1.15, whereas 
when 7 > 1.15 (and ih sufficiently small), it first decreases on the lower branch until it reaches a minimum and then 
increases along both branches. The magnitude of Q remains always small, however. For 7 = 0, Q approaches the 
finite limiting value Qn m ~ 0.124 on the upper branch, as wh tends to zero, independent of the isotropic horizon 
radius xh 6j. Apparently, this remains true when the dilaton is coupled. 

For comparison, we exhibit in Figs, la-c also the mass, the specific angular momentum, and the relative dilaton 
charge of embedded Abelian solutions with the same horizon radius (dotted curves) with Q = and P = 1. (Mass, 
angular momentum and dilaton charge of embedded Abelian solutions, possessing the same charge Q as the non- 
Abelian solutions and P = 1 are graphically indistinguishable from mass, angular momentum and dilaton charge of 
the Q = solutions, shown. Since the Abelian 7 = and 7 = V3 solutions are known analytically, their global 
charges along the upper branch are shown up to wh = 0.) As expected gg, mass, angular momentum per unit mass, 
and relative dilaton charge of the non-Abelian EYMD solutions, are close to mass, specific angular momentum, and 
dilaton charge of the embedded EMD solutions with Q = and P = 1. This is true, in particular, for solutions 
with small horizon radii ijj, where large deviations of these global properties from those of the corresponding Kerr 
solutions arise |6j]. The dilaton charge of the Abelian solutions, however, approach their limiting value D/j — on 
the upper branch with a common slope, different from the common slope of the non-Abelian solutions. 

The 7-dependence of the global charges is demonstrated further in Figs. 2a-b. Here we exhibit the mass M, the 
angular momentum per unit mass a = J/M, the relative dilaton charge -D/7, and the non-Abelian electric charge Q 
for non-Abelian black hole solutions with winding number n = 1, node number k = 1, horizon radius ih = 0.1, and 
wh = 0.020 on the lower branch (Fig. 2a) and on the upper branch (Fig. 2b) as functions of 7 for < 7 < 8. On 
the lower branch M, a = J/M, D/j, and Q decrease with increasing 7. In particular, we note that the non-Abelian 
electric charge Q passes zero at 7 ~ 1.5. On the upper branch only the relative dilaton charge D/7 decreases with 
increasing 7, whereas the mass M and the angular momentum per unit mass a increase with increasing 7. The 
non-Abelian electric charge Q exhibits a minimum on the upper branch. 

It is remarkable that the non-Abelian electric charge Q can change sign in the presence of the dilaton field 0. In 
contrast, in EYM theory the non-Abelian electric charge Q always has the same sign, depending on the direction of 
rotation [61I33T]. 

Thus we observe the new feature of rotating EYMD black holes that, for certain values of the dilaton coupling 
constant 7 and the parameters xh and wh, the non-Abelian electric charge Q of rotating black hole solutions with 
n = 1 and k = 1 can vanish. Cuts through the parameter space of solutions with vanishing Q are exhibited in Fig. 3a, 
where we show cjh as a function of xn for Q = solutions with dilaton coupling constants 7 = 1.5, V3, and 2. For 
fixed 7 these curves divide the parameter space into a region with solutions possessing negative Q and a region with 
solutions with positive Q. Negative Q solutions correspond to parameter values below the curves. 

The maximum of a curve coincides with w^ ax (:EH, 7). We note that, for fixed 7, the values of wh to the right of 
the maximum correspond to solutions on the lower branches whereas the values to the left correspond to solutions 
on the upper branches. In the limit 7 — > 7 m i n ~ 1-15 the curves degenerate to a point, implying that rotating Q = 
solutions exist only above 7 m i n - This is demonstrated in Fig. 3b, where the dependence of the maximal value of wh 
together with its location xh on 7 are shown, clearly exhibiting the minimal value of 7 for rotating Q — solutions. 

These Q = EYMD black holes represent the first black hole solutions, which carry non-trivial non-Abelian electric 
and magnetic fields and no non-Abelian charge. Perturbative studies [34| previously suggested the existence of rotating 
non-Abelian Q — black hole solutions in EYM theory, satisfying,however, a different set of boundary conditions. 
Such EYM black hole solutions could not be obtained numerically @, |2(j ■ 

As a consequence of their vanishing non-Abelian electric charge, these black hole solutions do not exhibit the generic 
asymptotic non- integer power fall-off of the non-Abelian gauge field functions For these solutions, the magnetic 
moment can be identified according to Eq. (|67fl . where fi x ~ C3 and fi z ~ — C%. The coefficient — C2 = |Cg| is exhibited 
in Fig. 3c for the sets of solutions shown in Fig. 3a. Note that the right end points of these curves correspond to the 
static limit. 

The corresponding curves for C3 are graphically almost indistinguishable. The difference between — C2 and C3 is 
due to the rotation, since in the static case C2 + C3 = 0. In Fig. 3d we show this difference for the same set of 
solutions. 

Having considered the rotating hairy black hole solutions for fixed horizon radius ih as a function of wh, we now 
keep wh fixed and vary the horizon radius. In Figs. 4a-d we show the mass M, the angular momentum per unit 
mass a = J/M, the relative dilaton charge D/j, and the non-Abelian electric charge Q of the non-Abelian black 
hole solutions as functions of the isotropic horizon radius xu for wh = 0.040 and 7 = 0, 0.5, 1, V3 and 3. (Solutions 
along the steeper branch are numerically difficult to obtain, therefore the D/j and Q curves are discontinued at 
in ~ 0.43.) For EYMD black hole solutions, as for EYM black hole solutions, there is a minimal value of the horizon 
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radius ^'"(wh), for a given value of wh- In particular, the limit x™ 1 — > is only reached for — > 0. We note that 
for KK black holes with Q = and P = 1 the family of solutions characterised by [w H , ^^"(wh)] tends to the static 
'extremal' solution as sh — ► 0. Anticipating that the non-Abelian black holes behave essentially analogously to the 
Abelian solutions, we expect to find the static regular EYMD solutions in this limit [||. We have no indication that 
rotating regular solutions could be reached in this limit |3J, |35j . 

Horizon Properties 

Let us next turn to the horizon properties of the rotating non-Abelian black holes. In Fig. 5a-d we show the area 
parameter x&, the surface gravity K sg , the deformation of the horizon as quantified by the ratio of equatorial to polar 
circumferences, L c /L p , and the Gaussian curvature at the poles K(0) as functions of w^, for black hole solutions with 
horizon radius xn — 0.1, winding number n = 1 and node number k = 1, and for dilaton coupling constants 7 = 0, 
0.5, 1, -\/3 and 3. Also shown are the horizon properties of the corresponding embedded Abelian black hole solutions 
with Q = and P = 1. 

The horizon size is quantified by the area parameter ita, shown in Fig. 5a. Starting from the value of the corre- 
sponding static non-Abelian black hole solution, the area parameter grows monotonically along both branches, and 
diverges along the upper branch. The horizon size of the corresponding Abelian black hole solutions is slightly larger, 
in particular along the lower branch. As illustrated, the difference decreases with increasing 7. The difference also 
decreases with increasing horizon radius ieh- 

The surface gravity n sg of the non-Abelian black holes is shown in Fig. 5b. Starting from the value of the corre- 
sponding static non-Abelian black hole solution on the lower branch in the limit ojh — ► 0, it decreases monotonically 
along both branches, and reaches zero on the upper branch in the limit ojh — ► 0, corresponding to the value assumed 
by extremal black hole solutions. The surface gravity of the corresponding Abelian black hole solutions is slightly 
smaller, in particular along the lower branch. The difference decreases with increasing horizon radius x h- W e recall, 
that for static non-Abelian black holes the surface gravity diverges in the limit xh — * independent of 7 36] , whereas 
for static Abelian black holes the limiting value depends on 7, and diverges only for 7 > 1, whereas it is finite for 
7 = 1, and zero for 7 < 1. 

The deformation of the horizon is revealed when measuring the circumference of the horizon along the equator, L e , 
and the circumference of the horizon along a great circle passing through the poles, L p , Eqs. I|79|l . The ratio L c /L p 
grows monotonically along both branches, as shown in Fig. 5c. As wh tends to zero on the upper branch, the ratio 
tends to the value L e /L p — n/(y/2E(l/2)) w 1.645, for both non-Abelian and Abelian black holes, independent of 7. 
Clearly, the matter fields lose their influence on the black hole properties, in the limit M — ► 00. On the lower branch 
the ratio L c / L p assumes the value one in the limit lur — > 0, corresponding to the value of a static spherically symmetric 
black hole. Overall, the deformation of the horizon of the non-Abelian black holes is close to the deformation of the 
horizon of the corresponding EMD black holes with Q = and P = 1. 

For rapidly rotating EYMD and EMD black holes, the Gaussian curvature of the horizon can become negative, as 
observed before for EM black holes In Fig. 5d we show the Gaussian curvature of the horizon at the pole K(Q) as 
a function of for non-Abelian and Abelian black holes. Starting with a positive curvature at the pole of the static 
solution, the curvature at the pole decreases with increasing angular momentum, crosses zero and becomes negative 
at some point along the upper branch. It then reaches a minimum value and starts to increase again, to become zero 
in the infinite mass limit, when cjh — * 0. 

In Figs. 6 we consider the deformation of the horizon of the non-Abelian black holes in more detail. We exhibit in 
Fig. 6a the angular dependence of the Gaussian curvature at the horizon for 7 = 1 black hole solutions with horizon 
radius sh = 0.1, winding number n = 1, node number k = 1, for several values of the angular momentum per unit 
mass a. When a > 0, the Gaussian curvature increases monotonically from the pole to the equator. For larger values 
of a, the Gaussian curvature is negative in the vicinity of the pole. 

Fig. 6b shows the shape of the horizon of the non-Abelian black holes, obtained from isometric embeddings of 
the horizon in Euclidean space [l6| . As pointed out in |16|. when the Gaussian curvature becomes negative, the 
embedding cannot be performed completely in Euclidean space with metric ds 2 — dx 2 + dy 2 + dz 2 , but the region 
with negative curvature must be embedded in pseudo-Euclidean space with metric ds 2 = dx 2 + dy 2 — dz 2 , represented 
by dashed lines in the figure. 

Limiting Abelian Solutions 

As seen above, the global charges and the horizon properties of the non-Abelian black hole solutions with n = 1 
and k = 1 are rather close to those of the embedded EMD solutions with Q = and P = 1. But also their metric 
and dilaton functions as well as their gauge field functions are impressively close to their EMD counterparts, except 
for those gauge field functions which do not vanish in the static limit. 
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With increasing node number k, the non-Abelian solutions get increasingly closer to the Abelian solutions, converg- 
ing pointwise in the limit k — *■ oo. In particular, the gauge field functions H 2 and H4 approach zero with increasing 
node number k in an increasing interval, while their boundary conditions #2(00) = 7/4(00) = (— l) fe enforce a value, 
differing from the value of the limiting Abelian solutions, H^(oo) = H^(oo) = 0. 

Convergence of the solutions with increasing node number k towards the corresponding EMD solutions is demon- 
strated in Figs. 7a-e for the metric function /, the dilaton function tj>, the gauge field functions H2, H3 and for the 
function xB>\ for non-Abelian black hole solutions with xh = 0.1, ojh = 0.020, n = 1, k = 1 — 3, and dilaton coupling 
constants 7 = and 1. From Fig. 7e we conclude that the limiting solution has vanishing electric charge. 

B. n > 1 Black Holes 

Here we give for the first time a detailed account of rotating n > 1 black hole solutions 7] . Whereas the rotating 
n = 1 black hole solutions reduce to spherically symmetric black hole solutions in the static limit, the rotating n > 1 
black hole solutions remain axially symmetric in the static limit 0, Q . 

Global charges 

Many features of the n > 1 EYMD black hole solutions, agree with those of the n = 1 EYMD solutions. In 
Figs. 8a-d we exhibit the mass M, the angular momentum per unit mass a = J/M, the relative dilaton charge -D/7, 
and the non-Abelian electric charge Q as functions of the parameter lor for black holes with horizon radius ih = 0.1, 
winding numbers n = 1 — 3, and node number k = 1 for the dilaton coupling constants 7 = 0, 1, and v3- Again, for 
7 = the relative dilaton charge D/j is extracted from the mass formula Eq. QJ, discussed below. 

The mass M shown in Fig. 8a and the angular momentum per unit mass a shown in Fig. 8b, increase monotonically 
and diverge on the upper branch in the limit ojh - * 0, as in the n = 1 case. Whereas the mass on the upper branch 
becomes (almost) independent both of the coupling constant 7 and the winding number n already for moderately 
large values for wh, the specific angular momentum becomes independent only of the coupling constant 7, but retains 
its dependence on the winding number n up to considerably smaller values of u>h- 

The relative dilaton charge shown in Fig. 8c, apparently approaches zero in the limit wh ™> 0, as in the case of the 
n = 1 black hole solutions. We note, however, that for winding number n — 3, the dilaton curve developes a minimum 
with a negative value of D, and only then bends upwards towards zero. The non-Abelian electric charge exhibited in 
Fig. 8d, increases with increasing n. As in the n = 1 case, Q apparently reaches a limiting value on the upper branch 
as wh — * 0. This limiting value increases with increasing winding number n and appears to be independent of the 
coupling constant 7. (For small values of wh on the upper branch, the numerical error becomes too large to extract 
the electric charge and thus the limiting values of Q.) For n — 3, the electric charge developes a maximum on the 
upper branch before tending to the limiting value. 

Comparing the global charges of these non-Abelian black hole solutions with those of the corresponding Q = and 
P = n Abelian black hole solutions, we observe that the discrepancies between Abelian and non-Abelian solutions 
become larger with increasing winding number n. As in the n = 1 case, the discrepancies decrease with increasing 
dilaton coupling constant 7, and with increasing horizon radius cch- 

Trying to find n > 1 black holes with vanishing non-Abelian electric charge, we constructed solutions for a wide 
range of parameters. However, in contrast to the n = 1 case, we did not observe a change of sign of Q for any set of 
solutions constructed. 

Horizon Properties 

Turning next to the horizon properties of the rotating non-Abelian n > 1 black holes, we exhibit in Figs. 9a-d the 
area parameter xa, the surface gravity K sg , the deformation of the horizon as quantified by the ratio of equatorial to 
polar circumferences, L e /L p , and the Gaussian curvature at the poles K(0) as functions of ojhj for black hole solutions 
with xh =0.1, n = 1 — 3 and k — 1, and for dilaton coupling constants 7 = 0, l,and y/3. 

As shown in Fig. 9a, the horizon size is monotonically increasing along both branches, and diverges as wh — > 
on the upper branch. Along the upper upper branch, the horizon size shows a distinct dependence on the winding 
number n, but little dependence on the coupling constant 7. The 7 dependence increases though with increasing n. 
The surface gravity decreases monotonically along both branches, and reaches a vanishing value on the upper branch 
as — > 0, as demonstrated in Fig. 9b. For black hole solutions with the same value of 7 but different values of n, 
we observe crossings of the lower branches. This is in contrast to the Abelian black hole solutions, where analogous 
crossings do not occur. It is thus a distinct feature of the non-Abelian nature of the black holes. 

The deformation parameter L e /L p is shown in Fig. 9c. For small values of wh along the upper branch, we observe, 
that the curves become (rather) independent of 7, but their slope keeps a distinct rt-dependence. All curves appear 
to have the same limiting value though, when — * 0, independent of 7 and n. (Unfortunately, the numerical error 



20 



for small wh along the upper branches is too large to extract the limiting value as wh — ► 0.) We expect this common 
limiting value to be L c /L p = 1.645, which is the value of the extremal Kerr black holes, Kerr-Newman black holes 
and Kaluza-Klein black holes. 

The Gaussian curvature of the horizon at the poles, exhibited in Fig. 9d, possesses negative values for fast rotating 
black holes, also for n > 1 black holes. The topology of the horizon of these n > 1 black holes is also that of a 
2-sphere, Eq. (JB2J. 

As expected, these properties of the non-Abelian black hole solutions are rather similar to those of their Abelian 
counterparts, though the differences become more pronounced as n increases. 

Limiting Abelian Solutions 

Considering the limit of node number k — ► oo and fixed winding number n as well as fixed 7, we observe that the 
black hole solutions with winding number n tend to the corresponding Abelian black hole solutions with vanishing 
electric charge Q — and with magnetic charge P — n. The convergence is observed for the global charges, and for 
the horizon properties, and pointwise convergence is seen for the metric, dilaton and gauge field functions. 

C. Mass Formula and Uniqueness 

The numerical stationary axially symmetric EYMD black hole solutions satisfy the mass formula, Eq. with an 
accuracy of 10~ 3 . So do the numerical EMD black hole solutions. EYM black holes are included in the limit 7^0, 
since D/7 remains finite. 

We demonstrate the mass formula for the non-Abelian black hole solutions in Figs. lOa-c, where we show the four 
terms in the mass formula, 2TS, 2fl J, -D/7, and 2ip c \Q, as well as their sum (denoted by M in the figures) as functions 
of the mass M. Black hole solutions with horizon radius xjj = 0.1, winding number n = 1, and node number k = 1, 
and dilaton coupling constant 7 = 1 are shown in Fig. 10a, n = 1 and k — 2 solutions in Fig. 10b, and n = 2 and 
k = 1 solutions in Fig. 10c. 

Reconsidering the uniqueness conjecture for EYMD black holes, one may try to obtain a new uniqueness conjecture 
by only replacing the magnetic charge P by the dilaton charge D or the relative dilaton charge -D/7. After all, such a 
replacement gave rise to the new non-Abelian mass formula, Eq. Q). However, in the case of the uniqueness conjecture, 
such a replacement is not sufficient. Indeed, black holes in SU(2) EYMD theory are not uniquely characterized by 
their mass M, their angular momentum J, their non-Abelian electric charge Q, and their dilaton charge D. This 
becomes evident already in the static case, where genuine EYMD black holes have J = Q = 0, and thus their only 
charges are the mass M and the dilaton charge D. In Fig. 11 we show the relative dilaton charge -D/7 as a function 
of the mass M, for the static non-Abelian black hole solutions with n = 1 — 3 and k ~ 1 — 3 for 7 = 1. Also shown 
are the corresponding limiting Abelian solutions with P — n. As seen in the figures, solutions with the same winding 
number n and different node number k do not intersect. However, solutions with different winding numbers n can 
intersect. The curve of n = 3, k = 1 solutions, for instance, intersects the curves of all n = 2, k > 2 solutions. 

To remedy this deficiency, we make use of the topological charge N, Eq. I|88|l . Recalling that for the non-Abelian 
black holes N — n, Eq. I]9ip. and setting N = for embedded Abelian black holes, we can attach to each black hole 
solution a topological charge. We then find, that all static SU(2) EYMD black hole solutions are uniquely determined 
by their mass M, their dilaton charge D (or the relative dilaton charge -D/7), and their topological charge N. We 
note, that for a given winding number N — n, the dilaton charge D (or the relative dilaton charge D/7) of the 
non-Abelian solutions increases monotonically with increasing fc, converging to the dilaton charge of the embedded 
Abelian solution with magnetic charge P = 11. 

This observation has led us to suggest a new uniqueness conjecture 0, stating that black holes in SU(2) EYMD 
theory are uniquely determined by their mass M , their angular momentum J, their non-Abelian electric charge Q, 
their dilaton charge D, and their topological charge N. Since the relative dilaton charge D/7 remains finite in the 
limit 7^0, this conjecture may formally be extended to the EYM case, by replacing the dilaton charge D by the 
relative dilaton charge -D/7. 

The uniqueness conjecture is supported by all our numerical data. To illustrate the validity of the conjecture also 
for stationary black holes, we show in Figs. 12a-c the relative dilaton charge -D/7 for several values of the specific 
angular momentum a for black holes with n — 1 — 3, k = 1 — 3 (7= 1). Also shown are the corresponding limiting 
Abelian solutions with P = n. For instance, we observe in Fig. 12a that the curve with n = 1, k = 2, and a = 0.25 
crosses the other k = 2 curves and also the k = 3 curves with a = and a = 0.10. However, since it does not cross 
the k = 3 curve with the same specific angular momentum a = 0.25, the conjecture is not violated. 
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V. 



CONCLUSIONS 



We have constructed non-perturbative rotating hairy black hole solutions of Einstein- Yang-Mills (-dilaton) theory 
and investigated their properties. The rotating black hole solutions emerge from the corresponding static hairy black 
hole solutions, when a small horizon angular velocity is imposed via the boundary conditions. The hairy black hole 
solutions are labeled by the node number k and the winding number n of their gauge field functions. We have extended 
our previous studies |6j, by considering black holes with higher winding number n, and by including a dilaton field with 
coupling constant 7. We have also provided a detailed analysis of the asymptotic expansion of black hole solutions 
with n = 1 and n = 3. 

The generic hairy black hole solutions possess a small non-Abelian electric charge, but vanishing magnetic charge. 
For fixed winding number n, the non-Abelian black hole solutions form sequences, labeled by the node number k. 
In the limit k — > 00, they tend to embedded Abelian black hole solutions with magnetic charge P — n and electric 
charge Q = 0. For 7 = the limiting solutions are Kerr-Newman solutions, for 7^0 they are rotating EMD solutions 
[tj flCX ■ A comprehensive study of the black hole solutions in EMD theory for general dilaton coupling constant 
will be presented elsewhere |17| |. 



The presence of the dilaton field has a number of significant consequences for the hairy black hole solutions. First, 
the presence of the dilaton field allows for solutions with vanishing electric charge Q. Thus there are rotating black 
hole solutions with no non-Abelian charge at all. These Q = solutions exist in a certain range of parameters for 
black holes with winding number n = 1, and node number k = 1. The asymptotic expansion of these Q = solutions 
involves only integer powers, in constrast to the generic non-integer power fall-off for solutions with Q ^ 0. We did 
not find corresponding Q = solutions with n > 1 or k > 1, when searching for them in a wide range of parameters. 

Second, the presence of the dilaton is crucial in deriving the mass formula |jj for rotating black holes, involving only 
global charges and horizon properties of the black holes. Here the derivation of the mass formula was presented in 
detail. The mass formula is similar to the Smarr formula for Abelian black holes, but instead of the magnetic charge 
P the dilaton charge D enters. This is crucial, since the hairy black holes have P = 0, but non-vanishing magnetic 
fields and thus a magnetic contribution to the mass, which is taken into account by the dilaton charge term. The mass 
formula also holds for EMD black holes. Involving the relative dilaton charge -D/7, the mass formula may formally 
be extended to the EYM and EM case, since the relative dilaton charge D/7 remains finite in the limit 7^0. 

Third, the presence of the dilaton allows for a new uniqueness conjecture for hairy black holes 0. Including the 
horizon topological charge N — n, our non-perturbative set of solutions supports the conjecture, that black holes 
in SU(2) EYMD theory are uniquely determined by their mass M , their angular momentum J, their non-Abelian 
electric charge Q, their dilaton charge D, and their topological charge N. Again, since the relative dilaton charge 
-D/7 remains finite in the limit 7^0, this conjecture may formally be extended to the EYM case, by replacing the 
dilaton charge D by the relative dilaton charge -D/7. 

We have not yet fully addressed a number of further issues for these rotating black holes. For instance, in the 
numerical part of the paper we have restricted the calculations to small values of the dilaton coupling constant 7. For 
larger values of 7, new phenomena might arise. We have neither considered in any detail the existence of stationary 
rotating 'extremal' EYMD solutions Q. Both these issues deserve further investigation. 

Most importantly, we have not addressed the existence of rotating regular solutions in EYMD theory. While 
perturbative studies indicated the existence of slowly rotating regular solutions |34j]. non-perturbative numerical 
studies did not support their existence [H,|2(j • Independent of the (non-)existence of slowly rotating solutions, there 
might however exist rapidly rotating branches of regular as well as black holes solutions, not connected to the static 
solutions. The investigation of such rapidly rotating solutions remains a major challenge. 



FNL has been supported in part by a FPI Predoctoral Scholarship from Ministerio de Education (Spain) and by 
DGICYT Project PB98-0772. 



The asymptotic expansion depends on two integers, labelling the black hole solutions, the winding number n and 
the node number k. We here explicitly show the dependence on fc, by choosing the boundary conditions according to 
Eq. (|36l) . Note, however, that all solutions can be brought to satisfy the same asymptotic boundary conditions by a 
gauge transformation (see the discussion in section II. C). Here we restrict to winding number n ~ 1 and n — 3. The 
analysis for n = 2 seems to be 'prohibitively complicated' [2(j because already in the static case the leading order 



Acknowledgments 



APPENDIX A: EXPANSION AT INFINITY 



terms involve non-analytic terms ~ log(x)/a; 2 |30| . 
n = 1 Solutions 

Generalizing our previous expansion for n = 1 [(| including the dilaton, we obtain 
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where a = a/9 - 4Q 2 and /3 = ^25- 4Q 2 . 



1 = solutions (n = 1, k = 1) 

Let us now consider the asymptotic expansion for the case Q = 0, and expand the constants Cj in powers of Q, 
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The Q-dependence of C4 and C5 is indeed observed numerically, as seen in Figs. 13a-b. For this figure we first 
obtained C20 and C30 from the Q — solution, yielding C 2 o = —25.014 and C30 = 25.019. The dashed lines correspond 
to the approximation, that only the first term in C4 and C5 is taken, using these values for C20 and C30. Then we 
extracted C4 and C5 from the function Hi for several values of 7 (equivalent to several values of Q, since and wh 
are fixed), where we employed C2 and C3, extracted from the functions H2 and H3 for the same values of 7. In the 
figure we observe good agreement between these points and the approximate theoretical curve. 

We also observe, that unlike the pure EYM case @, we may have C20 + C30 / 0, when the dilaton is included. 
Consequently, H3 can approach zero asymptotically with an integer power fall-off. Setting Cio = — 1/2A/, 2 , C\\ = 0, 
C20 = —b, C21 = 0, C30 = b, C31 = 0, and C%\ — 0, we obtain the relations for the EYM case, reported previously [||. 
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where e = \/49 - AQ 2 . 



APPENDIX B: EXPANSION AT THE HORIZON 
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Here we present the expansion of the functions of the stationary axially symmetric black hole solutions at the 
horizon xh in powers of 5. These expansions can be obtained from the regularity conditions imposed on the Einstein 
equations and the matter field equations: 
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APPENDIX C: ENERGY CONDITIONS AND SEGRE TYPE 



In this section we revisit the energy conditions for SU(2) EYMD theory, the results being valid for general gauge 
group SU(N). We also consider the Segre type of this type of matter. 

The dominant energy condition requires that "for every timelike vector V, the stress-energy tensor satisfies 
T^yV^V > and T^V V is a non-spacelike vector". Let us show that Eq. (J5J satisfies this condition. 

Let V be an arbitrary timelike vector. We may define its associated unit vector Eq via 



V = ^-V^E . (CI) 

By constructing an orthonormal basis with Eq as its timelike element, {Eq, E\, Ei, £"3}, one can reformulate the 
dominant energy condition as 

T( )(o) > |r (B)(6) | , Vo, b = 0,1, 2, 3, (C2) 

where the (•) index indicates the component in the orthonormal basis. 
In order to simplify the proof we separate Eq. JSJ into two parts, namely, 

T^=T™ + T™ , (C3) 



with 

ndU 



= d^d^ - ^g^d a ^d a <f = H^H V - ^g^H 2 , (C4) 
t ym = 2e 2 ^Tr(F MQ ^ 5 ^ - -g^F a0 F^) , (C5) 



where we have defined = i9 M $. 



A straightforward calculation performed in the orthonormal basis then yields 

T (o)(0) > \ T (a)(b)l T (0K0) > l T (aK&)l ' Vfl ' 6 = ' 1 ' 2 ' 3 ' ( C6 ) 

and from Eqs. I|C3() and l|C6|) we conclude that condition (|C2I) holds. 

Obviously, the weak energy condition is also satisfied, T(q\/q) > , Va, b — 0, 1, 2, 3. 

As for the strong energy condition, it can be formulated as follows: "A stress-energy tensor is said to satisfy the 
strong energy condition if it obeys the inequality 



T^VV > -TVpV , (C7) 



for any timelike vector V (T stands for the trace of the stress-energy tensor, T = T^)" . 

The Yang-Mills part of the stress-energy tensor, fulfils the strong energy condition trivially as it is a traceless 

tensor that satisfies the weak energy condition. The dilaton part, T°~, also satisfies the strong energy condition as 
can be easily shown in an orthonormal basis 

Note that T dil = T = -H 2 . Thus, Eq. (JC^ is satisfied by the full EYMD stress-energy tensor. 

Related to the energy conditions is the Segre type (or algebraic type) of the stress-energy tensor 18] . It consists of 
a study of the eigenvalue problem of that tensor, namely, 

T»„u u = \u» . (C9) 

It is well known that a general EM system has an algebraic type Al[(l 1) (1, 1)], with two double eigenvalues, 
opposite in sign. (We are not considering the null case). For an arbitrary EYM system the Segre type is that of a 
diagonalizable general type, Al[l 1 1, 1]. The non-Abelian nature breaks the degeneracy of the Abelian case. 

But the inclusion of the dilaton field also destroys that degeneracy Al[(l 1) (1, 1)] in the Abelian case, since the 
eigenvalues for a general EMD solution satisfy 

^1+^2 = — H 2 , 

A 3 + A 4 = , (CIO) 
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with H defined as before. It is no surprise then to find a generic type for the EYMD solutions. 
If we restrict ourselves to Eqs. (|15H and (|20|I - H22|I . we may write the eigenvalues in the form 



Ai 
A 2 
A 3 
A 4 



(T*H 


-T%) + 




- T"^) 2 + 47^7% 


(T*H 


-T%)~ 




- T^) 2 + 47^7% 


(T\- 


f T%)+, 




-T%) 2 +4T 8 r T r e 


(T r r - 


f T e e ) — a 




-T e e ) 2 +4T e r T I ' ( , 



(Cll) 



Note that this is just a consequence of the Lewis-Papapetrou form of the metric and the Einstein equations. By using 
an orthonormal basis it is easy to show that all the eigenvalues are real for EMYD solutions. 

The explicit expressions of these eigenvalues are rather complicated. To conclude this section we show the asymptotic 
behavior of the dimcnsionless version of these eigenvalues (n = 1) 



Ai 
A 2 
A 3 
A 4 



D 2 - Q 2 M(D 2 - 2Q 2 



2x 4 



1 



D 2 + Q 2 M(D 2 + 2Q 2 ) - 2 1 DQ 2 



2x* x 5 

D 2 -Q 2 M(D 2 - 2Q 2 ) 

2x~ A Hfi 

D 2 - Q 2 M(D 2 - 2Q 2 



2x A 



(C12) 



Obviously, we observe that for charged EYM solutions (D — 0), the algebraic type of the stress-energy tensor behaves 
asymptotically as that of a pure EM system. 
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Fig. la 
x H =0.1 , n=1 , k=1 
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FIG. la: The dimensionless mass M is shown as a function of luh for EYMD black holes with winding number n = 1, node 
number k — 1, horizon radius xh = 0.1, for dilaton coupling constant 7 = 0, 0.5, 1, \/3, and 3 (solid lines). The dimensionless 
mass of the corresponding EMD solutions with electric charge Q = and magnetic charge P = 1 is also shown (dotted lines) . 

Fig. lb 
x H =0.1, n=1 , k=1 
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FIG. lb: Same as Fig. la for the specific angular momentum a = J/M. 



30 



Fig. lc 

x H =0.1 , n=1 , k=1 
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FIG. lc: Same as Fig. la for the relative dilaton charge D/'y. 
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FIG. Id: The electric charge Q is shown as a function of uh for the same set of EYMD black holes shown in Fig. la. 



31 



Fig. 2a 
co H =0.020, x H =0.1 
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FIG. 2a: The global charges, M, a, D/'j, and Q, are shown as functions of 7 for EYMD black holes with winding number 
n = 1, node number k = 1, horizon radius :th = 0.1, and cjh = 0.020 on the lower branch (solid lines). The global charges of 
the corresponding EMD solutions with electric charge Q = and magnetic charge P = 1 are also shown (dotted lines). 

Fig. 2b 
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FIG. 2b: Same as Fig. 2a for black hole solutions on the upper branch. 
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Fig. 3a 

Q=0, n=1, k=1 

0.05 
0.04 
0.03 

x 

3 

0.02 



0.01 
0.00 

0.00 0.03 0.06 0.09 0.12 0.15 0.18 

X H 

FIG. 3a: Cuts through the parameter space of Q = black hole solutions with winding number n = 1, node number k = 1, 
and dilaton coupling constant 7 = 1.5, \/3, and 2. The curves are extrapolated to luh — xh = 0, as indicated by the dots. 

Fig. 3b 
n=1, k=1 

0.08 



0.06 

>t 0.04 
± 
S 

0.02 



0.00 

1.00 1.25 1.50 1.75 2.00 2.25 2.50 

y 

FIG. 3b: The value of the maximum of ujk and its location xk are shown as functions of 7 for the Q = EYMD black hole 
solutions with winding number n = 1 and node number k = 1. 
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Fig. 3c 
Q=0, n=1, k=1 



4.5 



3.5 



6 1 2.5 



1.5 



0.5 




0.00 0.03 0.06 0.09 0.12 0.15 0.18 



FIG. 3c: The coefficient — C2 = | C2 1 of the dimensionless magnetic moment fi is shown as a function of xh for the sets of 
EYMD back holes shown in Fig. 3a. 

Fig. 3d 
Q=0, n=1,k=1 
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FIG. 3d: Same as Fig. 3c for C 2 + C 3 . 
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Fig. 4a 
co H =0.040, n=1, k=1 
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FIG. 4a: The dimensionless mass M is shown as a function of xu for EYMD black holes with winding number n = 1, node 
number k = 1, wh = 0.040, for dilaton coupling constant 7 = 0, 0.5, 1, \/3, and 3. 

Fig. 4b 
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FIG. 4b: Same as Fig. 4a for the specific angular momentum a = J/M. 
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Fig. 4c 
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FIG. 4c: Same as Fig. 4a for the relative dilaton charge D/'y. 

Fig. 4d 
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FIG. 4d: Same as Fig. 4a for the electric charge Q. 
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Fig. 5a 
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FIG. 5a: The area parameter xa is shown as a function of u>u for EYMD black holes with winding number n = 1, node number 
k = 1, horizon radius xn = 0.1, for dilaton coupling constant 7 = 0, 0.5, 1, \/3, and 3 (solid lines). The area parameter of the 
corresponding EMD solutions with electric charge Q — and magnetic charge P = 1 is also shown (dotted lines) . 



Fig. 5b 
x H =0.1 , n=1 , k=1 
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FIG. 5b: Same as Fig. 5a for the surface gravity K sg . 
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Fig. 5c 
x H =0.1 , n=1 , k=1 
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FIG. 5c: Same as Fig. 5a for the ratio of circumferences L e /L p . 

Fig. 5d 
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FIG. 5d: Same as Fig. 5a for the Gaussian curvature at the poles K (0). 
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Fig. 6a 
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FIG. 6a: The Gaussian curvature at the horizon K(9) is shown as a function of the angle 6 for EYMD black holes with winding 
number n = 1, node number k = 1, horizon radius xu = 0.1, dilaton coupling constant 7 = 1, for the values of the angular 
momentum per unit mass a = 0, 0.257, 0.758, and 1.471. 



Fig. 6b 
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FIG. 6b: The shape of the horizon is shown for EYMD black holes with winding number n = 1, node number k = 1, horizon 
radius xu = 0.1, dilaton coupling constant 7 = 1, for the values of the angular momentum per unit mass a = 0, 0.76, and 1.47. 
Solid lines indicate embedding in Euclidean space, dashed lines in pseudo-Euclidean space. 
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Fig. 7a 
co H =0.020, x H =0.1, n=1 




FIG. 7a: The difference of the metric function / and the metric function f x of the limiting solutions is shown as a function of 
the dimensionless coordinate x for 8 = for the EYMD black hole solutions with 7 = (solid lines) and 7=1 (dotted lines), 
winding number n = 1, uih = 0.020 (lower branch), horizon radius xh = 0.1, and node numbers k = 1 — 3. 
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FIG. 7b: The same as Fig. 7a for the dilaton function <f>. 
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Fig. 7c 
co H =0.020, x H =0.1, n=1 
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FIG. 7c: The same as Fig. 7a for the function H-z- (iS~2oo = 0). 

Fig. 7d 
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FIG. 7d: The same as Fig. 7a for the function H z , but = 7r/4. 
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Fig. 8a 
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FIG. 8a: The dimensionless mass M is shown as a function of uh for EYMD black holes with n = 1 — 3, k — 1, = 0.1, and 
7 = 0, 1, and \/3. 

Fig. 8b 
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FIG. 8b: Same as Fig. 8a for the angular momentum per unit mass a. 
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FIG. 8c: Same as Fig. 8a for the relative charge D/7. 
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FIG. 8d: Same as Fig. 8a for the electric charge Q. 
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FIG. 9a: The area parameter xa is shown as a function of uh for EYMD black holes with n = 1 — 3, 
7 = 0, 1, and y/3. 



= 1, xh = 0.1, and 



Fig. 9b 
x H =0.1, k=1 



0.8 



0.6 



J 0.4 



0.2 



0.0 




0.000 0.009 0.018 0.027 0.036 0.045 

FIG. 9b: Same as Fig. 9a for the surface gravity K sg . 
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FIG. 9c: Same as Fig. 9a for the deformation parameter L e /L p . 
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FIG. 9d: Same as Fig. 9a for the Gaussian curvature at the pole K(0). 
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Fig. 10a 
x H =0.1 , n=1, k=1 , y=1 




FIG. 10a: Numerical check of the mass formula, Eq. ifT). for n = 1, k = 1, xh = 0.1, 7 = 1 black hole solutions. The terms in 
Eq. Q are plotted versus M, their sum coincides with M, within the numerical accuracy. 
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FIG. 10b: Same as Fig. 10a for n = 1 and fc = 2. 
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Fig. 10c 
x H =0.1, n=2, k=1 , y=1 




FIG. 10c: Same as Fig. 10a for n = 2 and k = 1. 

Fig. 11 
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FIG. 11: The relative dilaton charge D/7 is shown as a function of the mass for static non-Abelian solutions (7 = 1) and 
embedded Abelian solutions with Q = and P — n. 
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Fig. 12a 




FIG. 12a: The relative dilaton charge D/7 is shown as a function of the mass M for non-Abelian black holes with n = 1, 
k = 1 — 3 (7 = 1) and specific angular momentum a = 0.25, 0.1 and 0. Also shown are the embedded Abelian solutions with 
Q = and P = 1. 
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FIG. 12b: Same as Fig. 12a for n 



= 2 and P = 2. 
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FIG. 12c: Same as Fig. 12a for n = 3 and P = 3. 
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Fig. 13a 
co H =0.040, x H =1, n=1, k=1 
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FIG. 13a: The dependence of the asymptotic coefficient C4 on Q is shown for the black hole solutions with n = 1, k = 1, 
xh = 1, and luh = 0.040. The dashed line corresponds to the leading term in the expansion of C4. The dots correspond to 
numerical values of C4, extracted for several values of 7. 
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FIG. 13b: Same as Fig. 13a for C 5 . 



